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Abstract 

The general term of the Poincare normalizing series is explicitly con¬ 
structed for non-resonant systems of ODE’s in a large class of equations. 
In the resonant case, a non-local transformation is found, which exactly 
linearizes the ODE’s and whose series expansion always converges in a 
finite domain. Examples are treated. 

PACS code: 47.20.Ky. 


1 Normal forms and Poincare transformation. 

Let us consider a system of nonlinear ordinary differential equations: 

Xi = XiXifi{x), {i = l,...,N) (1) 

where the (purely nonlinear) functions fiix) are analytic in the real variables 
XI...XN. We suppose, for the simplicity of the purpose, that the linear part has 
been diagonalised. 

The study of the qualitative behavior of the system starts with a linear 
stability analysis of the fixed points |Q| . This determines the values of the control 
parameters for which the system bifurcates (i.e. for which some stable fixed 
point becomes unstable, or vice-versa.). A further analysis gives the nonlinear 
behavior of the solutions in the neighborhood of the hxed points. Among the 
methods to compute this last step, the most often used is the normal form 
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analysis 0, that we shortly describe here: consider a particular form of the 
system for which the origin is a fixed point; this implies the following form: 

ii = XiXi+'^ai{m)x"^, = (2) 

m 

where we used the multi-index notation, m = (mi, ...niN), ai{m) = ai{mi ,..., m^v) 
and x"* = The multiple sum is taken over the integers mi, 

such that X]^i —\ ^ 1^ 2, with rrii > 0. Consider now the near identity 
change of variables given by the following formal series: 

Xz = Vi + ^ h{m)y’^ (3) 

m,|m|>2 

This gives rise to a new system for the j/i: 

yt=Xzyi+ ^ Ci{m)y'^, (i = l,...,fV) (4) 

m,|m|>2 

Now, let us choose the bi{m) in such a way that this last system is as simple as 
possible. The ideal case would be Ci{m) = 0 for all i and all m, but it is easy 
to see that this will be possible only if no resonance condition is satisfied, i.e. if 
there exist no set of N positive integers with X]^i — 2, such that: 

N 

^rjAj-A*=0 (5) 

for some i. The reason is that, substituting (||) in (^) and requiring Ci{m) = 0, 
for all i and all m, we determine bi(rn) by a relation where the factor of bi(m) 
is precisely (X]j=i ~ ^ this quantity is equal to zero, bi{m) remains 

undetermined, and the system (^) will still be nonlinear, with: 

Ci(m) 7^ 0, if - Ai = 0 

= 0, otherwise (6) 

Equations with conditions (H) are called normal form. They form, in fact, 
a system equivalent to (H), where all the non-resonant monomial have been 
eliminated. The transformation is then called Poincare transformation, and 
the monomials remaining in (^, resonant monomials. 

In most cases, the general term of the series giving the Poincare trans¬ 
formation is not known, and, usually, it is computed order by order. When no 
resonance condition is satisfied, the solutions of (^) are simply exponentials of 
time, and (|^) provides a formal series of exponentials which is solution of (H). 
Moreover, the convergence of (||) is guaranteed in two cases [||: if the real part 
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of all the Xi are all of the same sign (Poincare domain). Or if the Xt belong 
to the Siegel domain, i.e. zero lies within the convex hull of Ai,Ajv in the 
complex A-plane, along with: 


N 


I- 

i=i 


for some c > 0, 

Nonetheless, there are domains of the complex A-plane for which the con¬ 
vergence of the Poincare series is not guaranteed. 

Even if the series is convergent, still one problem remains: the initial condi¬ 
tions of the new system have to be calculated in terms of the initial conditions 
of the old system i) and of its parameters Xi and ai{m). To investigate this, 
we have to invert the Poincare series, which is not an easy task. In practice, one 
inverts a truncated series. This leads to an error on the new initial conditions 
which may be excessively large and not controlled. 

In the resonant case, equation (^) is no longer linear, and contains in fact 
an infinity of monomials (since if one resonance condition is satisfied, so are an 
infinity of others). The use of the Poincare series to approximate the solution 
of (i is even more difficult, and suffers from the absence of general convergence 
theorems. 

From a geometrical point of view, the Poincare transformation defines a dif- 
feomorphic transformation to coordinates in the phase-space, in which the solu¬ 
tions are, in the non-resonant case, exponentials. The possible divergence in the 
non-resonant case, and the impossibility of eliminating the resonant monomials 
express the fact that the transformation leading to such a coordinate system is 
not always a diffeomorphism. 

The main goal of this paper is first to derive the general term of the Poincare 
series for an n-dimensional Lotka-Volterra (L.V.) system in the non-resonant 
case. That is, we find an exact solution to the recursion relation for the coeffi¬ 
cient of the Poincare series. In a second step, we build a Poincare-like transfor¬ 
mation such that: 

1- The system for the new variables is always linear, no matter if resonance 
conditions are satisfied. 

2- The initial conditions for the new system are known (we show that they 
may be chosen equal to the initial conditions for the original system.). 

3- The convergence of the series giving the transformation is always guaran¬ 
teed in a domain around the initial time to- 

The price to pay for such nice properties is to be found in the fact that 
this transformation is non-local as will be discussed later: it does not define a 
diffeomorphism. But still, it has many attractive features ranging from explicit 
computation of solutions to more theoretical questions in nonlinear dynamics. 

The restriction to the L.V. format may look rather restrictive. This is not 
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the case, since it is known that a large class of systems can be transformed into 
L.V. systems i i i i- We show in the next section how this can be done. 


2 The quasi-monomial formalism. 


We restrict our attention to nonlinear systems that can be written as: 

M N 

Xi= + {i = l,...N) (7) 

3 = 1 k=l 

No particular restriction is made neither on the Aij nor on the Bjk'. they may 
be real or complex numbers. M is the number of quasi-monomials 
appearing in (|^. In what follows, we suppose that M is finite. The results 
proposed in this work may be extended to the case of an infinite number M of 
quasi-monomials, provided some convergence conditions on the Aij are satisfied. 
The class (@) is thus quite general: many of the systems of interest in physics 
belong to this category. Systems of the type (|^) are called Quasi-Monomial 
(Q.M.) systems. 

The transformation to the L.V. format is made by adding M new variables 
to the system (Q). These are defined as follow: 

N 

XN+k = {k = (8) 

z=i 


Taking the derivative with respect to time, we find the {N -\- M) dimensional 
system: 


with: 


and: 

My 


N+M 

Xi — \iXi -f Xi ^ ^ ]\^lijXj (9) 

t=l 


Xi = ai, for 1 < i < N (10) 

N 


0, for 1 < z < + M; 1 < j 

A{j-n): for 1 < i < N] N ^ 1 < j < N ^ M (11) 

N 

^ ^ -^(z—w) 5 for A^ 1 ^ ^ ^ A^ H“ N H“ 1 ^ ^ A^ Af 
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The system (^) is of the form of the well known Lotka-Volterra equations, first 
introduced in theoretical ecology . 

Note that the linear term can be eliminated by adding one more variable 
which is put equal to one (xn+m+i = 0; XN+M+i{io) = !)■ We then have: 

N+M+l 
Xi — Xi ^ ^ 

i=i 

where: 

Mij = Mij, for l<i, j<M + N 

\i, ioi 1 < i < M + N] j = M + N + 1 

0, for i = M + TV + 1; 1 < j < M + iV + 1 

The main characteristic of the systems (||) and ( |l^ is that they present the 
lowest nonlinearity, i.e. the quadratic one. Moreover, the quadratic term is 
not the most general one (which would be of the form ^ NijkX^x^) but still 
permits complex behavior like chaos. 

It has been shown Q that the general structure for the Taylor series coeffi¬ 
cient can be obtained for the solution of One finds: 

W = 2^Ci{n) -1- 

n =0 

Ci(0) = Xi{t=tQ)=XiQ 

Ci(n) = XiO 

il ...in 

0^22 0 

We will show, in the next section, that the general structure of the Poincare 
series coefficient is very close to this one. In the last part of this paper, we also 
show that the Taylor series is a particular case of a class of series expansions of 
the solution which includes also the Poincare series. 


3 Structural analysis of the Poincare transfor¬ 
mation. 

We take as starting point (since the linear part plays an important role in the 
normal form analysis) the system: 

N 

Xi — \iXi ~h Xi ^ ^ JVlijXj , (z 1, A^) (1^) 
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Where, from now on, N stands for + M in the previous expressions. 

It is known, since the work of Carleman |^, that a nonlinear system 
can be viewed as an infinite-dimensional linear system. This can be realized by 
considering as new variables all the monomials one can build with products of 
positive integer powers of the Xi. Using the multi-index notation: 


Xm=X^=xrxr-xT ( 14 ) 

and taking the derivative of (|l^, one finds: 

N N N 

Xm = {J2 m + E(E 

k^l p^l 1^1 


where ep is a unit vector such that {ep)s = 5p^s (i.e. Xm+e^ = Xmi,...,mp+i,...,mN)- 
The infinite-dimensional linear system (H) is characterized by a triangular ma¬ 
trix R which is given by: 


N N N 

Rmp — E EE 

0=1 k=l 1=1 

For an original system &) that would be linear (that is, the matrix M 
vanishes) the system (El) would be purely diagonal. This implies that the 
Poincare transformation on (H, for non vanishing matrix M, but in absence 
of any resonance, corresponds to the diagonalisation of the infinite-dimensional 
matrix defined by the system (HI)- 
Consider now the relation: 


L 


mp 


= 5 


m,p 


+ 


Rmk ( 1 - ^mk)Lkp 

Rpp Rnim 


(16) 


where, once again, the indices are multiple (the sum over fc is a multi-sum over 
the N indices ki, ...,k]si running from 0 to oo, and 6m,p stands for Smi,pi---SmN ,pn)- 
We claim that since Rmp is triangular (that is: Rmp = 0 if there exists at 
least one integer k between 1 and N, such that nik > Pk) and if Rmm ^ Rpp for 
all m,p with m p, then {Lmp) diagonalises {Rmp)- 

The condition Rmm ^ Rpp for m ^ p implies that the relation 


N 

^kTk = 0 

k^l 

is satisfied only if = 0 for all k and restricts thus the system dl) to the non¬ 
case. That {Lmp) diagonalises {Rmp) means here that, considering the inverse 
operator defined by: 


E ^mk^kp — y] L 


mk^ 


- -1 


= <5. 


mpt 


k 


k 


(17) 
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we have: 


^ ^ ^mk^^oLop — Rrap^n 


(18) 


The question of the existence of the inverse operator L~^ is obvious. This 
operator represents the inverse of the Poincare transformation. The latter is a 
diffeomorphism. Hence, when it exists, so does its inverse. 

Let us now prove this proposition: 


Multiplying both sides of by {Rpp — Rmm + Sm,p) one finds: 


k 


R^nkLkp — LmpRpp ^mpi^Lmp 1 ) ^mp {Rpp — R mm) 


(19) 


We will now show that the last two terms of ( 0 ) are equal to zero. It is clearly 
the case for Smp{Rpp — Rmm)- For Smp{Lmp — 1) we now demonstrate that, 
thanks to the fact that Rmp is triangular, so is also L^p, and all the elements 
Lmm are equal to one. To show this, consider the series obtained by iterating 
This will be a series of powers of Rmpi^ — ^mp) (with coefficients depending 
on m and p). Now, the elements of the power of i?mp(l ~ 5mp) for which 
p < m + k are all equal to zero. This implies that the only contribution to 
Lmm comes from the first term of the series, that is Smp- So, Lmm = 1, and 
5mp{Lmp — 1) = 0. Multiplying then (^^ by L~^ and summing over m, we find 
the announced result. 

Taking this result back to the original L.V. system, we can build the Poincare 
series. We insert in (p^, the matrix R defined by (jl^): 

and, writing 0) for m = Ci, we find after some simple algebra: 


Xi 


oo N N N 

—l^Ti )] ['^*71 — 1 ) * ■ ■ ('^*1 • ■ ■ H“ )] 2 /*! yZ 2 ’ 2/^71 


( 20 ) 


where the term corresponding to n = 0 is, by convention, set equal to one. 

This result can be readily extended to the system by using the defi¬ 
nitions 0) and (|Tl| ). If the system (i is polynomial, the obtained series is 
the Poincare series. But the class of systems given by ( 0 ) also contains non¬ 
polynomial systems. In these cases, the series which is derived is more general 
than the Poincare transformation of ®- 

Comparing this to the Taylor series, we see that the two structures are very 
close. This is not surprising: the Poincare series is, in fact, the Taylor series in 
which an infinity of resummations have been performed. 

Accordingly with what we said about the Poincare series, (^^ is not always 
convergent, as we shall see in the examples treated later. Moreover, this result 
is not valid in the resonant case, but in some sense, can be extended to it 


as 


we now see. 
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4 The Poincare transformation revisited. 

Let us define the N variables: 



where the Xi satisfies Taking the derivative of © with respect to time, 


we see that the differential system for Ui is always linear, iii = AiWi, no matter 
if resonance conditions are satisfied, and independently of the value of 7 . As 


announced, the definition (|^) is non-local. This is expressed by the presence 
of the integral: the value of Ui{t) depends on the values of Xi{t') for t' between 
7 and t. The inverse transformation (p^) reveals this nonlocality in an even 


clearer way. 

Choosing 7 equal to the initial time to, Ui{t) and Xi{t) have the same initial 
conditions. Moreover, ([2l|) is a recursion relation defining Xi{t) in term of Ui(t). 


It is possible to prove that generates the series: 


00 N N N 


X. = “.wEEE-i; 





[Ui^ (ai)Ui2 (On)] 


( 22 ) 


Here again, the explicit expressions of (^) and (|^) for the system are 
obtained using ( p^ an d (pi] ). 

Comparing this to (i^ we see that the \i appearing in the denominator, 
disappeared in the integrals of (©). Moreover these integrals are finite, since 
the integrands are exponentials, and the domains of integration are finite. Using 
the explicit form of Ui(t), we have: 


00 N N 




n=02l=l = l 


n=02l=l = l 



The integrals appearing in (^^ can be performed in closed form by using Carl¬ 
son special functions |^, but we do not need this result for our present purpose. 


Moreover, it can be proven that this series always converges on a finite domain 
around to, even in the resonant case. These results will be detailed in a forth¬ 



coming publication . 
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can then be performed order by order, and we see that Ui(t) will be equal to yi(t) 
if 7 is fixed in such a way that the constant term appearing after performing 
the integral (i.e. the primitive of evaluated at t = 7 ), is equal to 

zero. If the real parts of the A^’s are all positive (resp. negative) we have to 
put 7 = —00 (resp. + 00 ), since the constant term is a sum of exponentials of 
7 . But if there exist A^’s whose real parts are of opposite sign, there will be 
exponentials with positive arguments, and others with negative arguments. In 
this case, if 7 tends to infinity (no matter if it is +00 or — 00 ), some exponentials 
will be divergent. Identically, in the resonant case, some monomials built on 
the variables ?/i’s are constant, and no value of 7 will fulfill the requirement. 

In fact, the correspondence between (^) and (p^ ) is more subtle. A bijective 
correspondence can be established between the terms appearing in ( p 0 |) and 
rooted tree graphs. In the graph formalism, ( 1 ) represents a relation between 
these graphs |^. 

An easy way to see why ( |^) can be divergent even in the non-resonant case, 
and in the same blow, why (| 2 l]) is always convergent in a certain domain, is 
the following: series (p0| ) is a Taylor series expansion (in powers of yi) of the 
Poincare transformation around the origin in the y-space. But nothing secures 
that the trajectory of the solution will go through the origin. On the other 
hand, by definition, (0 guarantees that the series reduces to its order zero at 
t = to- In other words, the solution in the j/-space is, at t = to, precisely at the 
point around which the series is built. 

Finally, as announced, when all the A^ tend to zero, the integrals in ( p^ 
generate the factors {t — to)^/n\, and ( p^ is precisely the Taylor series given 
above, whereas series ( p0| ) for A^ tending to zero diverges. 

Thus, although the transformation (EiD is non-local, its structure explains 
many aspects related to the existence and convergence of the Poincare trans¬ 
formations. Moreover, it leads to useful algorithms for explicit calculation of 
solutions. 


5 Example: Ricatti projective systems. 

As an example, we now treat a particular case of the system (^^: a class of 
V-dimensional, integrable systems, called projective Ricatti systems Hi 0: 



The usual Poincare series (0) gives: 
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Thanks to the fact that the general solution of the system is known, we can 
determine the initial conditions for the yi. This leads to: 


Xi{t) 




i-E 


N 


Cj! XjQ 
^3 


OO 


E( 

n=0 


CjXjo \,(t-tn) 
1 _ CkXko 

^ l^k=l Afc 


)n 


It is clear that, for some choice of the Ci, Xi, and Xio, this series will be divergent 
for any t. For example, in two dimensions, if we choose Ai real and positive 
(Ai = Z), A 2 = —Ai, XiQ = 3:20 = S; and Ci = 1 = —C 2 , the series becomes: 


Xi{t) 




1 - 2s/l 


\ 10^/7 ^ 


n—Q 


1 - 2s// 


And the radius of convergence is given by the inequality: 

I ^ ch{l{t-to)) \< 1 
l/s — 2 

If I > 4s, there exist no t for which the above relation is satisfied. 
Consider now the series (^3|): in this case, it gives: 


CXD N 

n—0 j—1 


A 3 ^ 


1 ))” 


which is always convergent on a domain ]to — T,to + t[. Moreover, if one of the 
Ai’s tends to zero, the limit can be taken without difficulty (and is finite), which 
is not the case for the usual Poincare series given above. 
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